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Abstract. We report the detection of anomalous brightness fluctuations in the multiple image Q0957 + 561 A,B 
gravitational lens system, and consider whether such anomalies have a plausible interpretation within the frame- 
work of cosmic string theory. We study a simple model of gravitational lensing by an asymmetric rotating string. 
An explicit form of the lens equation is obtained and approximate relations for magnification are derived. We 
show that such a model with typical parameters of the GUT string can quantitatively reproduce the observed 
pattern of brightness fluctuations. On the other hand explanation involving a binary star system as an alternative 
cause requires an unacceptably large massive object at a small distance. We also discuss possible observational 
manifestations of cosmic strings within our lens model. 
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1. Introduction 

Recent observations of the Q0957+561 A,B gravitational 
lens system show unexpected synchronous (without the 
expected time delay) fluctuations of brightness of the two 
quasar images. An ordinary binary star system, which 
might theoretically explain such fluctuations, would be 
too massive and close to us, i.e., would be clearly visi- 
ble, which is not the case. Therefore, we attempt to ex- 
plain these data by lensing on a cosmic string, particu- 
larly on a loop of string. The existence of cosmic strings is 
predicted by particle physics (Vilenkin & Shellard 1994) 
and gravitational lensing effects are a promising signa- 
ture of these astrophysical objects. Sazhin et al. (2003) 
claimed the detection of the first case of cosmic string lens- 
ing. We demonstrate here another possible signature of a 
string: microlensing by oscillating loops of cosmic strings, 
which results in quasiperiodic fluctuations of the observed 
brightness of the source. In section [21 we discuss the ob- 
servational data, and in section [3] a quantitative model of 
string lensing is elaborated. In section0]an explanation of 
the observational data is presented, with discussion and 
conclusions in section [5] 



2. Observations of an anomaly in the Q0957+561 
A,B gravitational lens system 

The Q0957 system was the first discovered multiple im- 
age gravitational lens system, and already at the time of 
its discovery in 1979 (Walsh et al. 1979) it was under- 
stood that it was extremely important to astrophysics. 
Measurement of the time delay between fluctuations in 
the two known images would allow determination of the 
Hubble constant from simple theory, independent of un- 
certainty in local distance estimates for the supernovae 
and Cepheid variable stars (Refsdal 1964) 



Thus from the 

time of discovery, monitoring of the brightness of the two 
images, separated by 6 arcsec, was undertaken so that 
the quasar's intrinsic brightness fluctuations could be rec- 
ognized in the two images separately, and a time delay 
measured. 
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With the Schild and Cholfin (1986) measurement of 
time delay (including numerous subsequent refinements; 
see Colley et al (2003) for a summary) it was soon rec- 
ognized that there were differences between the time de- 
lay corrected brightness curves, although the basic pattern 
could be easily recognized. The differences were attributed 
to microlensing by individual massive objects, presum- 
ably stars, in the lens galaxy (Schild & Smith ITI)91|) . The 
prospect that such microlensing might allow detection of 
any baryonic missing mass objects justified intensive mon- 
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Fig. 1. Brightness of the two quasar images displayed with 
no correction for gravitational lens time delay. The bright- 
ness of quasar image A (upper record with square data 
markers), has been fitted with a sine curve having 0.04 
magnitude amplitude. The lower record, with triangular 
markers, appears to have the same sinusoidal brightness 
curve with zero lag, even though at most epochs data for 
the gravitationally lensed images show a lag of 417 days. 



itoring campaigns, and in the 24 years since discovery the 
source has been consistently observed on more than 1500 
nights. 

Such monitoring reveals two principal components in 
the quasar's brightness fluctuations: a component due to 
intrinsic quasar brightness fluctuations, first seen in im- 
age A and then seen 417.1 days later in image B, and a 
microlensing component arising in only one image compo- 
nent due to individual stars along the A or B image line 
of sight. 

We now illustrate what appears to be a third compo- 
nent of quasar brightness fluctuations, seen in the com- 
bination of Figs, n 121 an d El Here we plot the measured 
brightnesses of the two quasar images as measured during 
the 1994-1995 and 1995-1996 seasons. In Fig. Q] no cor- 
rection has been made for time delay; the plotted bright- 
ness measurements in magnitudes are shown for the Julian 
dates of observation. The plotted symbols are the size of 
the typical 1 — <r error bars previously established for this 
data set (.006 and .007 mag for images A and B). In Fig.^ 



Fig. 2. Quasar brightness displayed for measured time de- 
lay. The upper record shows the same data and fit for 
image A as displayed in Fig. ^ The lower data mark- 
ers (triangles) are the brightness measurements for image 
B measured 417 days (the gravitational lens time delay) 
later, but with 417 subtracted from the Julian dates for 
plotting. If the image A brightness fluctuations are intrin- 
sic to the quasar, they should be seen also in image B 417 
days later, but the two are seen not to match as well as 
the lag comparison in Fig. ^ 



a sine curve has been fitted to the A image data but not 
B image to allow the eye to judge whether there appears 
to be a repeating pattern of fluctuations for lag. 

Fig. ^ shows the unexpected result that a short- 
duration oscillation of 4% amplitude and periodicity of 
approximately T ss 100 days was seen for approximately 
400 days. The amplitude of these fluctuations is well above 
the known errors of the measurements. The error esti- 
mates originally attributed to these data by Schild l|1995f) 
have been confirmed from subsequent analysis by Colley 
and Schild ( 1999 ), and the entire data set has subsequently 
been re-reduced by Ovaldsen et al. I|2()()3 ). who also call 
attention to the observed correlation for lag. 

The correlation for lag is not perfect, as would be ex- 
pected, since several processes are simultaneously causing 
brightness fluctuations. Microlensing can impose a ran- 
dom pattern of fluctuations with durations ranging from 
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Fig. 3. Quasar brightness with 417-day lag for the oppo- 
site image pair. The lower record is the Fig. ^ data for 
image B, with the data for image A measured 417 days 
previously. The agreement is seen to be poor, especially 
near the end of the observational period, even though the 
theory of gravitational lenses shows that the time delay 
must produce agreement for 417 days. If the measured si- 
nusoidal oscillation seen in both images A and B (Fig.[TJ is 
a chance coincidence of two quasar oscillations separated 
by 417 days, there must be agreement for both A data 
with corrected B data (Fig. |2J) and B data with corrected 
A data (this figure). Comparing Fig. ^ with Figs. |21 and [3] 
shows best agreement for lag, contrary to gravitational 
lens time delay theory. 



1 day to decades. An example of an event with 0.01 mag 
amplitude and 12 hour duration has been given by Colley 
and Schild ( 2003 ) . A wavelet analysis of the long bright- 
ness record by Schild (1999) shows that events on time 
scales of 1 and 60 days have typical amplitudes of 0.01 and 
0.08 magnitudes, respectively. Yet the fluctuation pattern 
is sometimes seen to effectively stop, as reported in Colley 
et al. PKEty . 

If the observed fluctuations are caused by random mi- 
crolensing events, it would be unexpected for them to be 
apparently in phase. We have not yet devised a statisti- 
cal test defining the significance level or error limits on 
simultaneity because the basic statistical process is non- 
Gaussian and has not yet been simulated. Moreover, any 



statistical analysis cannot be perfect in the presence of 
the usual stochastic microlensing variation taking into ac- 
count the limited time interval where the anomalous effect 
has been observed. If the fluctuations are intrinsic to the 
quasar, and seen simultaneously by some highly improba- 
ble coincidence, they must be seen in the observations of 
the preceding and following years, as illustrated in Figs. [3 
and|3 Thus we show in Fig. |2] that if the B data of Fig.^ 
are compared to A image data measured 417 days pre- 
viously, the fluctuations are probably not seen. Similarly 
we illustrate in Fig. [3] that if the A image pattern from 
Fig-Dis compared to B image 417 days later, the pattern 
is again not seen. If the brightness fluctuations arc intrin- 
sic to the quasar and seen simultaneously as in Fig. ^ by 
chance, they must also be seen in the lagged data for the 
opposite quasar image; thus they would be seen in both 
Fig. 0and Fig. 

The importance of these observations relates to the 
fact that there should be no causal connection between 
brightness fluctuations seen simultaneously. If the fluctua- 
tions were due to the quasar's intrinsic brightness changes, 
then they should be seen at the measured time delay, 
which is 417.09 ± 0.07 days (Colley et al. I271TT3T) . If they 
were produced in proximity to the lens galaxy at redshift 
0.37 (the quasar redshift is 1.41) they should similarly be 
seen with a large time delay. Only fluctuations produced 
locally (i.e., between the lens galaxy and the observer but 
close to the observer) can be observed to be simultaneous. 

This problem is even more serious because of the rela- 
tively large separation of the two quasar images, 6.2 arc- 
sec on the sky. Supposition that the above oscillations 
are induced by orbiting of binary stars leads to anoma- 
lously large masses of the components, as shown in subsec- 
tion ^^ Therefore we consider below the possibility that 
the oscillations are due to time variations of the gravita- 
tional field of a cosmic string. 



3. Gravitational lensing by a cosmic string 

3.1. Cosmic string characteristics 

Cosmic strings are linear defects that could be formed at a 
symmetry breaking phase transition in the early Universe 
(Vilenkin & Shellard lTTJQll) . 



A horizon-sized volume at any cosmological time t 
should contain a few long strings stretching across the 
volume as well as a large number of small closed loops. At 
the moment of creation t, typical loop length is I ~ act, 
i.e., about a of horizon size ct. During the string evolution, 
loops constitute some fixed part of total string network; 
this scaling results in the following loop number density 



ni(t) 



\ct) 



-3 



(1) 



Typically a is determined by the gravitational back- 
reaction, so that a ~ k g G[j,/c 2 — k g e, where k g ~ 50 
is a numerical coefficient, G is gravitational constant, c is 
the speed of light, fx ~ rj 2 /hc 3 is the mass per unit length 
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of string, Ti is the Planck constant, and r\ is the symme- 
try breaking scale of strings. For a Grand Unified Theory 
(GUT) string ? ?GUT ~ 10 25 eV and p G ^ ~ 10 22 g/cm, 
e GUT r^ 10~ 6 and «gut ~ 10 . 

The loops oscillate and lose their energy, mostly by 
gravitational radiation. For a loop of length I, the oscil- 
lation period is Ti = I /2c and the lifetime determined by 
gravitational radiation losses is 77 ~ l/c(k g Gp/c 2 ). 

Gravitational lensing by cosmic strings has been con- 
sidered by many authors (see references in Vilcnkin 
& Shellard (fTTTMjl and de Laix & Vachaspati l|T35S|l ). 
Straight cosmic strings have a distinctive feature: they 
produce two identical images. However they cannot ex- 
plain the oscillatory character of our data. Therefore we 
consider gravitational action of cosmic string loops, which 
cause effects similar to ordinary oscillating systems (bi- 
nary stars and others), but are more massive and move 
with relativistic speeds. 

3.2. Lensing by oscillating loops 

De Laix & Vachaspati l|1996[) considered in detail, lens- 
ing by cosmic string loops. Here we use their approach 
for the interpretation of the observed oscillations. In the 
simplest idealized case of a circular loop, with oscillations 
reduced to variations of loop radius, they find that the 
image brightness of a point source will not oscillate if the 
loop does not overlap the source. Consequently we should 
take an asymmetric loop to explain the observed oscilla- 
tions. We consider a maximally asymmetrical string con- 
figuration in the form of a rotating double line segment of 
length 2R lying transverse to the line of sight and having 
coordinates: 



i?cos (ct/R) sin ct, 
R sin (ct/R) sin a, 
0. 



(2) 



This is a particular solution from the family of known solu- 
tions to string equations (Turok 1984, Vilenkin & Shellard 
I1994[) . Here axis X3 is directed to observer, xf tr , x s 2 tr are 
coordinates of the loop in the lens plane, a indicates po- 
sition on the string and varies from to 2-7T and t is the 
time. This configuration is a limiting case of asymmetric 
loop; in a more realistic case a loop will have an ellipsoidal 
form with large eccentricity. 

The lens equations can be obtained from general result 
of de Laix & Vachaspati ijl99fi|l . After some calculations 
taking into account the particular solution @ we have: 

j7-yi =Xi- q s R 2 sgn(xi) x 
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(4) 



where q s = 8^-7^ J ^, 3 ^' (D s , Di and D[ s are distances 
from us to source plane, to lens plane and from source to 
lens plane respectively), x\ = x\ cos (ct/R)+x 2 sin (ct/R), 
x 2 = x-i cos (ct/R) — x\ sin (ct/R) (where X\, x 2 are co- 
ordinates in the lens plane) and y\ = y± cos (ct/R) + 
y 2 sin (ct/R), y 2 = y 2 cos (ct/R) - y\ sin (ct/R) (where y\, 
y 2 are coordinates in the source plane). It may be shown 
that all the relations under the root signs are non-negative. 
Magnification of a point-like source by such a string is 



l- q 2 R\x 2 



9 ct 

-AR (xi cos V x 2 sin 

R 



(R 2 



-4? 



Ct 2 \3/2 

R/ 



(5) 



If q 8 ~ 1 and R y = RD S /Di <C p = (y\ + y 2 ) 1/2 , we can 
solve approximately the lens equations (l.'il H and obtain 
the magnification 



+ 



= 1 + 
3g 2 i?6 



„2 7?4 
Qs U y 



p 6 



((v? 



y 2 )cos — 



2yxy 2 sin 



2ct\ 
~RJ' 



(6) 



This yields amplitude of source brightness fluctuations as 
follows 



A? 



6g*R« _ ZMtt 2 G 2 p 2 6\ 



(f 



C^j 



(7) 



where 9j = p/D s is the angular impact distance of the 
line-of-sight with respect to the center of the loop, 9r = 
R/Di is half of the visible angular size of the loop. 

3.3. Lensing by a binary system 

Now we consider for comparison the gravitational lensing 
by a binary system of two equal point masses M orbiting 
their center of mass with the period T. Further, r is half 
of the distance between the masses and us = n/T. The 
magnification of a point-like source by such lens system is 
(Schneider, Ehlers & Falco UMfy : 



-1 2 4 ( 1 2 1 2 1 i\i 
n,= 1 - q h r ( (xf + x 2 + r ) - 

-Ar 2 (xi smut - x 2 cos cut) 2 J / ((x\ + x\ 
-4r (x\ cos uit + x 2 sm cot)' 



„2\2 



(8) 



where q\> — 8 ^ M D j^ D i and the other values are defined as 
in the previous section. 
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Analogously to the previous case we obtain approxi- 
mate formulae for magnification 



m 



1 



9 A 



4ggr« 



6g2 r 6 



((2/1 ~ v\ ) cos 2wi + 2yiy 2 sin 2uit) 



(9) 



(where r y = rD s / D{) and for the amplitude of fluctuations 
in case of circular motion of the binary system: 



Am : 



12ggrg _ 1.2 x 10 6 efPfDJ 
p e c^T^D 2 



(10) 



where 0i — p/D s and r = r/Di (angles in radians). Here 
we have used the relation GM = 4ir 2 r 3 /T 2 for the circular 
motion. 



4. Explanation of the experimental data 

Finally we apply the above calculations to explain the ob- 
served brightness oscillations. These oscillations are nearly 
sinusoidal, their period is approximately 100 days and 
their amplitude is about 4% of the quasar image bright- 
ness. At least three oscillations were observed during the 
period of the observations. We consider the possibility that 
this phenomenon is caused by the cosmic string loop pass- 
ing through the neighborhood of images A and B at a 
small distance from the observer. Obviously to fit the ob- 
servational data described in the section [2] we are forced 
to restrict the parameters of our model. Also, because the 
synchronous oscillations have been observed within a lim- 
ited time interval, we include into the consideration the 
motion of the loop. At that the number of observed os- 
cillations (3-4) restricts a transverse velocity of the loop 
to the values wi,W2 < 0.1c, but leaves a considerable 
freedom for the velocity component V3 along the line of 
sight. In this case the only correction for the lens equa- 
tions, as can be shown, is to change the parameter q s by 
q^qsil + vs/c)- 1 . 





-4. -2. 0. 
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Fig. 4. Oscillations of quasar image brightness caused by 
the loop depending: a) on mass per unit length p, in g/cm, 
b) on transverse velocity of loop vj_ in units of the speed 
of light, c) on minimal distance po between image and 
center of the loop in units of R, d) on angle (po between 
loop direction and direction from the loop center to the 
image at t = in degrees. 



4.1. String-induced and binary star-induced oscillations 
of quasar brightness 

As we mentioned above, the period of observed oscillations 
T k, 100 days is related to the string length I — 2ttR as 
R = cTi/ir. Taking into account relativistic motion of the 
string along the line of sight we have TJ = T(l — t^/c 2 ) 1 / 2 . 
Equation can be rewritten as: 



Am 



■»■«( ^ r7 "' ; " 



/' 



1.5" 

2 



10 22 g/cm 



1 



V3 



(11) 



To provide almost equal amplitude of the brightness 
variations in both images the loop should fly close to the 
mid-point. In this case 9j « 3", i.e. about half of the 



distance between the images A and B. For numerical esti- 
mates we put, e. g., v « 0.7c. In order to have 3 oscillations 
with possible phase shift (see Fig. EJ we need v\ « 0.03 
and v 2 «0.11 (Fig. El. 

In order to have quasi-sinusoidal variations, Or must 
be considerably smaller than the angular distance between 
images A,B of the Q0957+561 (otherwise there will be 
sharp spikes and/or discontinuities in the dependence of 
brightness upon time) and Or cannot be too small lead- 
ing to large loop mass in virtue of equation Ijlll) (to avoid 
a large monopole input into the effective lens potential 
leading to additional amplification and corresponding un- 
observed slow increase and decrease of image brightness 
superimposed on smaller oscillations due to loop rota- 
tion). Therefore we should take Or « 1.5" and conse- 
quently Di « 3 kpc. From equation (|ll|l we can find that 
p, « 4 x 10 20 g/cm for observed amplitude Am w 0.04. 
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Fig. 6. a) View of the lens galaxy microlensed by the string loop at t — 0. A and B show the position of quasar images. 
Loop positions in different moments of time are indicated as well. Simultaneously b) the caustics (shown sideways) 
and the boundaries of the image doubling zone (upper and lower) in the plane of the galaxy and c) the critical curves 
around the loop edges in the lens plane are presented. 




3. T 



Fig. 5. Oscillations of quasar image magnification pre- 
dicted by the cosmic string model. Upper and lower curves 
are shifted up and down by 0.05 and fitted to image A and 
B brightness records, respectively. 



More accurate values of loop parameters follow from nu- 
merical solution of equations © - © without assumption 
y/yf + V% ^ -^y For explanation of observations we need 
<7 S = 1.3 and consequently [i ~ 8 x 10 21 g/cm (see Fig.©. 
Remarkably, this value is close to one predicted by the 

GUT, [i ~ /i G UT- 

In the alternative case of lensing by a binary system 
we obtain from equation l|10|) : 



0.04 



T 



100 days/ V 1 - 5 " 

6/ ,4 



1.2 pc 



(12) 



In order to explain the observed fluctuations the require- 
ment on the trajectory should be the same as in the 
loop case: 8i « 3", i.e. about half of the distance be- 
tween images A and B. The minimum distance from us to 
this system must be 1.2 pc, the orbital radius should be 
6 r Di « 1.8 a.u. and the masses of the components should 
equal 78Mq to supply 4% amplitude fluctuations of pe- 
riod 100 days. At such distance the stars would be observ- 
able. For larger distances to the binary star the masses of 
the components need to be larger as well. Therefore we 
consider such a binary star model to be unacceptable. 

4.2. Influence of the loop on the brightness of the 
lensing galaxy 

Let us now consider the influence of a loop on the bright- 
ness of the lens galaxy which is visible between the images 
of the quasar and makes a [ 3 %, 15 % } contribution to the 
observed brightness in the [ A, B ] apertures, respectively 
(Colley & Schild EH Fig. 3). For a galaxy at distance 
approximately 1" from image B with Gaussian brightness 
distribution (a = 2"), we obtain from numerical calcula- 
tion that relative oscillations of galaxy brightness in the 
measured apertures are equal to 9% and 6%. This corre- 
sponds to 1.6% and 1.1% of total signal in measured aper- 
tures A and B, which include quasar images. Oscillations 
will be superimposed on a background of monopol com- 
ponent with amplitudes of 12% and 14% of galaxy bright- 
ness. This corresponds to 2.2% and 2.5% of total bright- 
ness variations. A view of the lens galaxy as modified by 
the loop at t = is shown in Fig. 0a. 

The loop overlapping a galaxy can result in significant 
microlensing effect. Caustics and critical curves are shown 
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in Fig.|S|(b and c). When a star of radius i? s is located on 
the caustic near its edge, the star's brightness is increased 
by a factor of 3 x 10 5 (R s /Rq)~ 1 / 2 according to equation 
for magnification near a straight caustic (see Schneider, 
Ehlers, & Falco 1992]). The corresponding increase of the 
galaxy brightness will be about 3 x 10~ 7 i.e., unobserved 
in our case. 

The largest magnification is expected for stars cross- 
ing the cusp of the caustic. Our calculations show that 
magnification of a star with R s — Rq, is equal to 
m w 1.1 x 10 9 (0.1% relative to the galaxy bright- 
ness). The visible speed of the cusp motion through 
the galaxy plane is v cusp sa 2 x 10~ 3 pc/s. In 
one pixel of a telescope image plane, light is col- 
lected from about (D g 9 pixc i) 2 n stal ~ 5 x 10 8 (6> pixc i/0.1") 2 
stars (?i s tar ~ 1000 pc~ 2 is the surface density of stars, 
D g ps 1.4 Gpc is the distance to the galaxy). Therefore 
during about 10- 8 (6>p ixe i/0.1")~ 3 w 4 x 10~ n s the bright- 
ness of the star is larger than the brightness of the pixel 
without lensing. 

For telescope integration time Tj nt , the av- 
erage magnification of the star in the cusp is: 
< m >= 1.2 x 10 4 (r in t/ls)~ 2/3 w 400. The typical 
distance between projections of stars in the galaxy plane 
equals ^ ss « 0.03 pc. Therefore such flashes will be 
repeated approximately every l ss /v, 



cusp 



10 seconds. 

In our case the galaxy brightness is about 20% of the 
quasar B image brightness and therefore the maximum of 
these flashes will be approximately 2 x 10~ 4 of the image 
brightness, which is below observational precision. So 
lensing of the galaxy is observable only from distortion of 
the total galaxy image, as shown in Fig. a. 



5. Conclusions 

Our main motivation to apply the cosmic string model for 
explanation of brightness oscillations in question is their 
specific characteristics. It is difficult to propose a less ex- 
otic model, such as a double star model to explain the 
observed oscillations. We have shown (see subsection 14. l|l 
that in case of a double star the masses of the components 
must be of order of 78Mq at 1.2 pc distance in order to 
meet all observational requirements. For larger distances 
the masses of the components should be larger as well. 
This seems to be unacceptable. 

On the other hand the property of fast oscillations 
is typical for cosmic strings. To show viability of the 
string interpretation of observed oscillations we have cho- 
sen a particular "degenerate" (highly asymmetrical) ana- 
lytical solution of string equations. This solution is a limit- 
ing case of strongly elongated rotating loop configuration 
with sufficient quadrupole moment. Of course, more real- 
istic case should include additional modes of loop oscilla- 
tions in order to avoid self-intersections and annihilation. 
Nevertheless, even in more general cases solutions with 
sufficient quadrupole moment can provide the same level 
of brigthness variations. 



The results presented here show that loops of cosmic 
strings supply quantitative explanations of synchronous 
variations in the two images of the gravitationally lensed 
quasar Q0957+561 A,B. The derived value of the string 
parameter fi ~ 10 22 g/cm lies just in the theoretically 
preferable range. Atypical (with small probability of real- 
ization, but not impossible) in our model are the distances 
to the loops and their sizes - both are about 10~ 3 — 10~ 4 
of statistically expected values. The reason for this is the 
relatively short observed period, only of order 100 days, 
which limits the length of the loop. Moreover, the ob- 
served angular separation of the images predicts the dis- 
tance to the loop, while we fixed its size according to the 
period of flux fluctuation. Consequently in another hy- 
pothetical object with different observational parameters, 
more typical loops will probably work. Of course, a sin- 
gle event need not follow statistical rules, especially since 
our observational sampling rate might strongly favor this 
particular specimen. Moreover, some physical mechanism 
might cause the concentration of loops in galactic halos. 
Therefore, further observational efforts towards uncover- 
ing more objects with similar properties are extremely im- 
portant. 

Searches for brightness oscillations, similar to those de- 
scribed above, can be a promising way of discovering the 
gravitational signatures of cosmic strings. 
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